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ABSTRACT 

We present the first detection of molecular gas cooling CO emission lines from ordinary massive galaxies at 
Z = 1.5. Two sources were observed with the IRAM Plateau de Bure Interferometer, selected to lie in the 
mass-star formation rate correlation at their redshift, thus being representative of massive high-z galaxies. Both 
sources were detected with high confidence, yielding L' co ~2x 10 10 K km s _1 pc 2 . For one of the sources we 
find evidence for velocity shear, implying CO sizes of ~ 10 kpc. With an infrared luminosity of LpiR ~ 10 12 L«, 
these disk-like galaxies are borderline ULIRGs but with star formation efficiency similar to that of local spirals, 
and an order of magnitude lower than that in submm galaxies. This suggests a CO to total gas conversion factor 
similar to local spirals, gas consumption timescales approaching 1 Gyr or longer and molecular gas masses 
reaching ~ 10 n MQ, comparable to or larger than the estimated stellar masses. These results support a major 
role of in situ gas consumption over cosmological timescales and with relatively low star formation efficiency, 
analogous to that of local spiral disks, for the formation of today's most massive galaxies and their central black 
holes. Given the high space density of similar galaxies, ~ 10~ 4 Mpc~ 3 , this implies a widespread presence 
of gas rich galaxies in the early Universe, many of which might be within reach of detailed investigations of 
current and planned facilities. 

Subject headings: galaxies: evolution — galaxies: formation — cosmology: observations — galaxies: star- 
bursts — galaxies: high-redshift 



1. introduction 

Galaxies in the distant Universe were forming stars with 
much higher rates than today (Madau et al 1996; Elbaz et al 
2002; Le Floc'h et al. 2005). For galaxies of a fixed stellar 
mass, e.g. of the order of M* ~ 10 10_11 Mq, the typical star 
formation rate (SFR) has been steadily increasing with respect 
to local values, up to ~ 30 times higher at z = 2 (Daddi et al. 
2007a; D07a hereafter). This increase can be accounted for ei- 
ther by postulating a higher gas content or a higher star forma- 
tion efficiency, or both. Galaxy mergers and interactions, ex- 
pected to be much more common in the distant Universe than 
today, can increase the density of gas in galaxies, and thus 
help transform gas into stars more rapidly and with a higher 
star formation efficiency. The relative importance for massive 
galaxy formation and assembly of the hierarchical merging of 
small sub galactic units, versus the in situ formation of stars 
through gas consumption is largely unknown observationally. 

The molecular gas content of galaxies is best probed 
through the detection of redshifted cooling emission lines of 
CO. Combining observations for local and distant galaxies, 
all the way from spirals to (ultra) luminous infrared galax- 
ies (LIRGs, Lfir £ 1O U L , and ULIRGs, L Fm > 1O 12 L ), a 
nonlinear correlation is observed (e.g., Solomon & van den 
Bout 2005) between CO luminosity (proportional to molecu- 
lar gas content) and FIR luminosity (proportional to the SFR), 
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with more luminous sources having higher FIR to CO lumi- 
nosities, i.e. higher efficiencies than spiral disks in converting 
gas into stars, likely due to their higher gas densities (Tac- 
coni et al. 2006; Downes & Solomon 1998). Mapping the 
physics of molecular gas in representative galaxies at high 
redshifts would be a major step forward into understanding 
the processes regulating galaxy formation, and is a key target 
of future astrophysical facilities such as ALMA. 

However, if the local correlations hold at high redshifts, cur- 
rent techniques and sensitivities allow us to detect CO lines, 
and thus study the molecular gas content, only for extreme ob- 
jects such as bright (Ssso^,,, > 5 mJy) submm-selected galax- 
ies (Genzel et al. 2004; Neri et al. 2003; Greve et al. 2005). 
(SMGs), quasars (Walter et al. 2003; Maiolino et al. 2005), 
or in the case of strong gravitational lensing. Indirect evi- 
dence has been mounting very recently that this might not be 
the general case. Vigorous starbursts at high redshift appear 
to be the ordinary star formation mode in massive galaxies. 
The existence of tight correlations between SFRs and galaxy 
masses (Elbaz et al. 2007; Noeske et al. 2007; D07a) and 
the predominance of ULIRGs inside mass limited samples of 
M > 5 x 10 10 M@ galaxies (Daddi et al 2005; D07a; Caputi et 
al 2006) suggests that vigorous starbursts were typically long 
lasting among ordinary massive galaxies in the distant Uni- 
verse, which therefore might be hosting large amounts of gas. 

2. co observations of BzK galaxies 

Motivated by this evidence, we observed two representa- 
tive massive galaxies at high redshift with the IRAM Plateau 
de Bure Interferometer, searching for the redshifted emission 
lines from the CO [2-1] transition at rest frame 230.539 GHz, 
with the aim of constraining their molecular gas content. The 
two galaxies, BzK-A\ll (J123626.53+620835.3, VLA radio 
frame) and fiz^-21000 (J123710.60+622234.6), at redshift 
z = 1 .465 and z = 1 .522, respectively (measured from Keck 
DEIMOS spectroscopy), were selected in the Great Observa- 
tories Origins Deep Survey Northern field (GOODS-N) us- 
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FIG. 1. — Top panels: CO [2-1] emission line detections in two ordinary, 
massive high-z galaxies. We show (uncleaned) velocity averaged spatial maps 
of CO emission, from -350 km s _I to425kms _I (BzK - 4171, left) and be- 
tween -575 km s -1 and -25 km s _1 and from -950 km s and -800 km s 
(BzK — 21000, right). The contours are from -3 to +6<r, in steps of la 
(0.1 mJy/beam). The field of view is 1'. The size of the beam is shown in the 
lower-left. Bottom panels: position-velocity maps, extracted along the major 
axis direction measured on the HST+ACS UV rest frame imaging. The spa- 
tial field of view is ~ 9". Red/white colors are for positive pixels. Contours 
are at the level of 1 mJy/beam. 
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FIG. 2. — We show the CO spectra in velocity space, with respect to the 
zero velocity offsets (dashed vertical lines) defined from the optical redshift 
measured from the [Oil] emission lines (we assume an intensity weighted 
wavelength of 3727. 3A for the doublet). Best fitting Gaussian profiles are 
shown for the emission of the 2 targets. 



ing the BzK technique (Daddi et al 2004b), which provides 
complete ^-limited samples of normal, massive galaxies at 
1.4 < z < 2.5. SFRs of about 100-150 M yr~ r were de- 
rived by averaging estimates based on the VLA radio detec- 
tions of 40 ± 5/iJy and 49 ± 7/iJy (using the radio-LpiR cor- 
relation, Condon et al. 1992) and UV estimates (correcting 
for dust reddening using the Calzetti et al. 2000 law), which 
agree within 0.2 dex or better. For BzK — 4171 the 24^m 
estimate (based on the Chary and Elbaz 2001 luminosity de- 
pendent templates) also agrees within this factor, while for 
BzK - 21000 the 24^m estimate of SFR is 2-2.5 larger than 
the radio and UV ones, mild evidence for possible mid-IR 
excess usually linked to obscured AGN activity (Daddi et al. 
2007b; D07b hereafter). Stellar masses of 0.5-0.8 x 10 n M e 
were derived from optical and UV rest frame photometry, fol- 
lowing the method of Daddi et al. (2004b). SFRs and stellar 
masses are given here for the Chabrier (2003) stellar initial 
mass function. The two target galaxies lie in the mass-SFR 
correlation established at their redshifts (D07a), and are thus 
ordinary massive high-z sources with ongoing star formation. 

The sources were observed using newly refurbished re- 
ceivers for 5.7 and 7.5 hours during April 24-25th 2007 for 
BzK -411 1 and Bz^-21000, respectively, during stable and 
high quality conditions. A compact configuration-D with 6 
antennas was used. Additional 7.5 hours of poorer quality 
weather data in configuration-D and with 5 antennas only 
were obtained for BzK-4111 on August 1st 2007, reduc- 
ing the noise by only 25%. Standard data reduction was 
performed using the GILDAS software. Both sources ob- 
served were detected in CO with high confidence (Fig. [T] 
top panels), with a significance of 6.8cr and 8.9c. Fluxes 
are 0.60 ±0.09 Jy km s" 1 and 0.85 ±0.10 Jy km s _1 , for 
BzK-4111 and BzK-21000, respectively. Fig. |2] shows the CO 



spectra of the two BzK galaxies. The CO lines wavelengths 
are in agreement with the Keck redshifts, within 50 km s _1 
or better. The CO positions agree with the VLA positions 
within sa 0.5". BzK-4111 is best fitted by a double horn pro- 
file, indicative of a rotating disk and consistent with the UV 
morphology of a nearly edge on disk. BzK -21000 shows a 
narrower velocity profile with only a hint of double-horns. 
Fig- [U (bottom panels) shows position-velocity diagram, de- 
rived using a pseudo-slit oriented along the major axis de- 
fined in the HST+ACS imaging. BzK-4111 is unresolved, 
while BzK-21000 shows evidence for shear, overall consistent 
with rotation in a disk. As a cross-check, we separately fitted 
velocity channels bluer and redder than the systemic veloci- 
ties. For BzK-4111 the positions of the blue and red horns are 
separated by 0.5 ±0.5". For Bz^-21000 we find 1.3 ±0.3" 
separation (11 ±2.5 kpc; 4.3a), definying a position angle of 
— 10 ± 17 deg, in very good agreement with the UV position 
angle, confirming that the CO emission is extended. 

3. GAS RICH MASSIVE GALAXIES AT z ~ 1.5 

These results unveil a new and different galaxy formation 
mode for the most luminous sources like LIRGs and ULIRGs. 
In Fig. [3] we compare CO to FIR luminosities of our BzK 
galaxies to local and distant sources. We adapted published 
luminosities for local (Solomon et al. 1997; Yao et al 2003) 
and distant sources (Greve et al. 2005) to WMAP cosmologi- 
cal parameters (Spergel et al. 2007). We assume here roughly 
constant brightness temperature (see Braine & Combes 1992) 
and thus conversion ratios among the different CO tracers to 
be of order unity, i.e. CO[l-0]-CO[2-l]~CO[3-2] (WeiB et 
al 2005). Compared to local ULIRGs, these ordinary mas- 
sive z — 2 galaxies have similar bolometric luminosities but 
larger CO luminosities by a factor of 3-5. On the other hand, 
SMGs have estimated Lfir ^ 1O 13 L0 and therefore (Kenni- 
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FIG. 3. — Comparison of molecular CO luminosity (L' c0 ) to far-infrared luminosity (Lpm, integrated between 40 and 500^m). Quantities for the 2 BzK 
galaxies presented in this work are shown as large filled circles. The left panel shows the comparison of the 2 luminosities. The dotted line is the best fit relation 
to spirals (triangles), LIRG/ULIRGs (crosses) and distant SMGs (empty circles). The solid line is for a linear L' c0 to Lfm correlation normalized to the ratio of 
local spirals. The local LIRG with high L' c0 content, behaving similarly to distant BzKs, is VII Zw 31 (plotted as an asterisk; Downes & Solomon 1998), which 



of redshift. The dashed horizontal lines show the semi-interquartile range of local spirals. 



cutt et al. 1998) SFR > 1000M© yr~\ one order of magni- 
tude higher than our two targets, but the observed CO lumi- 
nosities are comparable. The two BzK galaxies are outliers 
of the Lfir versus Leo correlation defined by local spirals, 
LIRGs/ULIRGs and distant SMGs (Fig. |3). Instead, they lie 
close to a linear Lfir to Leo relation normalized to local spi- 
rals. The observed ratio in the two massive BzK galaxies is 
an order of magnitude lower than that of SMGs at similar red- 
shifts, implying much lower star formation efficiency. This is 
consistent with the fact that SMGs are outliers of the mass- 
SFR correlation (D07a), probably due to the higher efficiency 
in forming stars, for a similar stellar mass and CO luminosity. 
The Lfir to Leo ratios of the two BzK galaxies are within the 
semi-interquartile range of values from local spirals (Fig. [3] 
right panel). No object with similar Lfir to Leo ratio and large 
Lfir was previously known in the distant Universe (Fig. |3). 
This implies that these distant BzK galaxies behave as scaled 
up versions of local spirals, with proportionally higher SFR 
and molecular gas content for a fixed galaxy stellar mass. 

We use a conversion rate (Kennicutt et al. 1998) Lfir/L = 
0.8 x \0 w SFR [M s yr _1 ], adapted to a Chabrier (2003) stel- 
lar initial mass function and in the assumption, appropriate 
for these galaxies, that Lfir is about 80% of the bolometric 
luminosity. In general, we can write L co = M^/Xco, where 
Mh2 is the total molecular gas mass (including He) and the 
conversion factor (Downes & Solomon 1998; Solomon & van 
den Bout 2005) is X co = 4.6 M (Kkm s _1 pc 2 )" 1 for the 
Milky Way and as low as X C o = 0.8 M (K km s~ ! pc 2 ) -1 
for ULIRGs and likely for SMGs. The high detected CO lumi- 
nosities imply that giant molecular gas reservoirs are present, 
of order of 2 x 10 M , for a ULIRG-like X C o ~ 1, compa- 
rable to SMGs and larger than in local ULIRGs. The typical 
gas consumption timescales (Pig. [3j are about 2-A x 10 7 yr for 
SMGs, a factor of w 2 larger for local ULIRGs and reach 2- 
3 x 10 8 yr for the massive BzK galaxies at z = 1 .5, in the case 
of Xco ~ 1 ■ The SFR could last for even longer than that, 



as further gas accretion will be substantial during such long 
timescales. Hydrodynamical simulations from Keres et al. 
(2005) predict gas accretion rates of about 10-70 M yr~' for 
similarly massive galaxies at 1 .4 < z < 2.5 (i.e., for dark mat- 
ter halo masses of ~ 1O 12 ~ 13 M ; Hayashi et al 2007). This is 
the same order of magnitude of the estimated SFR in massive 
z ~ 2 galaxies. These are likely lower limits to the SFR dura- 
tion and gas masses as we are entirely neglecting atomic gas. 
As the star formation efficiencies are consistent with those of 
local spirals, in analogy to the latter it might be that Xqo > 1 
or even that the Milky Way conversion factor applies to these 
massive BzK sources. 

This is further supported by the evidence that, unlike local 
ULIRGs (Downes & Solomon 1998) and distant SMGs (Tac- 
coni et al 2006), these two BzK galaxies, like z > 1 .4 massive 
star forming galaxies (Daddi et al 2004a; Bouche et al 2007) 
in general, are not compact merger-induced starbursts. From 
detailed Sersic fitting of their UV rest frame light profiles 
from Hubble Space Telescope imaging data obtained with Ad- 
vanced Camera for Surveys (Giavalisco et al. 2004), we de- 
rive half-light radii of 4.4 and 5.5 kpc (proper size, or about 
0.5", accurate to better than 5%). Both galaxies have disk- 
like (exponential) morphology with Sersic index n = 1.3 and 
0.6 (accurate to better than 20%). BzK-4171 appears indeed 
as a nearly edge on disk, consistent with the Sersic fitting. 
BzK-21000 has a clumpy structure, still consistent with a disk 
but also with the possible presence of minor merging events. 
The UV rest frame sizes are likely reasonably accurate mea- 
surements of the starbursts sizes, because for these sources 
reliable SFRs measurements can be obtained from the UV 
(Daddi et al 2005; D07a; Dannerbauer et al 2006; Reddy et 
al 2005), i.e. their UV emission is not entirely opaque, while 
local ULIRGs and distant SMGs are optically thick in the UV 
(Champan et al 2005; Goldader et al 2002). This is confirmed 
by the evidence for extended CO emission, at least for BzK- 
21000. The inferred sizes are a factor of 2-3 larger than dis- 
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tant SMGs (Tacconi et al 2006) and 10 times larger than lo- 
cal ULIRGs (Downes & Solomon 1998). This converts into 
1-2 order of magnitudes lower gas densities and justifies the 
relatively modest star formation efficiencies. Independent es- 
timates of the molecular gas masses and thus SFR timescales, 
can also be obtained. First, if we apply the Schmidt-Kennicutt 
law (using its latest derivation from Bouche et al 2007, that 
includes also distant SMGs), using the UV size for the typ- 
ical gas and SFR size, we estimate M gas ~ 10 M© for each 
galaxy. Second, from the CO velocity widths and (UV and/or 
CO) size measurements of the galaxies we estimate dynam- 
ical masses of about 1—1.5x10 Mq, within the half light 
radii, including gas, stars and dark matter. Given the es- 
timated total stellar masses of ~ 0.5-0.8x10 M Q , roughly 
half of which presumably lying within the half-light radii, it 
looks plausible that w lO n M of gas can be present. We con- 
clude that it would take up to 1 Gyr to consume all of the gas 
in these typical, massive high-redshift galaxies, assuming the 
star-formation rate remained roughly constant. Clearly, the 
star formation episodes could last longer since additional gas 
is accreted over such long timescales. 

4. DISCUSSION 

In order to evaluate the implications of these results for 
massive galaxy formation, recall that the two observed galax- 
ies are representative of massive galaxies at high redshifts. 
They were selected for CO observations simply due to the 
availability of known spectroscopic redshifts. The two galax- 
ies lie on the SFR to stellar mass correlation inferred at 1 .4 < 
Z < 2.5, and their colors and morphologies are fully typical of 
massive z ~ 2 galaxies. Star forming galaxies in that redshift 
range with similar stellar masses are common in the distant 
Universe, with space densities of order of 10~ 4 Mpc~ 3 , quite 
a bit larger than the 2 - 6 x 10~ 6 Mpc~ 3 observed for SMGs 
by Chapman et al. (2003). Clearly, these results point to the 
existence of a much larger population of gas rich galaxies in 
the distant Universe than previously thought. 

We notice that the far-IR luminosity estimates for our ob- 
jects are based on the UV, radio and mid-IR luminosities, and 
are solid given that self-consistent measurements of SFR ac- 
tivity have been shown to be obtainable for these kind of ob- 
jects at high redshift (D07a). It is very unlikely that we are 
strongly underestimating their star formation activity (which 



would lead us to overestimate the duration timescales). In- 
stead, the Lfir for SMGs from Greve et al. (2005) is based 
only on the 850/im flux and is likely more uncertain. Us- 
ing radio and Spitzer, Pope et al. (2006) suggest lower LpiR 
for SMGs. This could imply that SMGs are closer in their 
SFR modes to massive BzK galaxies, but would change little 
about our conclusion because it is the latter which represent 
the dominating star formation mode at high redshift. 

The large molecular gas fraction in these galaxies, possi- 
bly up to 50% or more, would likely give rise to large scale 
instabilities, with the formation of clumps (e.g., Immeli et 
al. 2004; Bournaud, Elmegreen & Elmegreen 2007). This 
could explain the clumpy structure of these two galaxies, and 
in general of a large fraction of the high-z massive galaxies 
(e.g., Daddi et al. 2004a), and the high velocity dispersions 
observed in many of the high-z massive disk galaxies (e.g., 
Genzel et al. 2006; Forster-Schreiber et al. 2006). 

The long duration timescale for star formation activity im- 
plies, as a corollary, long duration activity for the galactic nu- 
clei. In fact, it has been recently unveiled by D07b that up 
to 20-50% of massive galaxies in 1.4 < z < 2.5 likely con- 
tain strongly obscured but relatively luminous AGNs. Some 
of the large gas reservoir is likely slowly funneled, in some 
way, into the central black hole with similar timescales to star 
formation, at a rate roughly consistent with what is required 
to assemble the correlation between stellar and black hole 
masses as we observe today (D07b). Our results underline 
the major role of gas consumption over long timescales and 
with low efficiencies, as opposed to rapid and strong merger 
driven bursts, as a major growth mode for both stellar mass 
and black holes in the distant Universe. 
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